Precise Measurement of the Pion Axial Form Factor in the 7r+ — > 6+1/7 Decay 
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We have studied radiative pion decays vr"*" ~* e'^vy in three broad kinematic regions using the 
PIBETA detector and a stopped pion beam. Based on Dalitz distributions of 42,209 events we 
have evaluated absolute tt evy branching ratios in the three regions. Minimum fits to the 
integral and differential {E^+ , Ej) distributions result in the axial-to- vector weak form factor ratio of 
7 = Fa/Fv = 0.443(15), or Fa = 0.0115(4) with Fv = 0.0259. However, deviations from Standard 
Model predictions in the high- i?^, /low- kinematic region indicate the need for further theoretical 
and experimental work. 

PACS numbers: ll.30.Rd, 11.40.-q, 13.20.Cz, 14.40.Aq 



In the Standard Model description of radiative pion de- 
cay TT"*" — > e^v^, where 7 is a real or virtual photon (e~^e~ 
pair), the decay amplitude M depends on the vector V 
and axial vector A weak hadronic currents . Both cur- 
rents contribute to the structure-dependent terms SDy 
and SDa associated with virtual hadronic states, while 
only the axial-vector current contributes to the inner 
bremsstrahlung process IB. Thus, it is convenient to 
write the decay amplitude as a sum: M = Mib + Msd- 
The IB amplitude is 2, 3J 



MiB = - I — -j= — Uruee'^ t 
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where p, k and q are the pion, electron and photon four- 
momenta, respectively, e and nie are the electron charge 
and mass, Gf is the Fermi coupling constant, Vud is the 
CKM quark mixing matrix element, while is the pion 
decay constant. 

The structure-dependent amplitude is parameterized 
by the vector and axial vector form factors, Fy and Fa'- 

Msd = j=-e e^^{l - 75)1^ 

X [Fve^iuaTP^'q' +iFA{g^ii,pq~PuqiL)\ ■ (2) 

The conserved vector current (CVC) hypothesis 
relates Fy to the 7r° lifetime yielding Fy = 0.0259(5) 
which agrees with the relativistic quark model and chiral 
erturbation theory . Chiral symmetry calculations 0, 
' 1 yield Fa in the range 0.010-0.014. 
The combined tt ev^ event count of all previously 
published experiments is less than 1,200 events, while 



the overall uncertainties of the parameter 7 = F^ 
extracted from data range from 12% to 56% |10llll 
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In this Letter we present a first analysis of the 7r+ — > 
e+z^7 (7re27) events recorded with the PIBETA detector 
in the course of a new measurement of the 7r+ — > Tr^e^i/ 
(tt/?) branching ratio J^ilJ from 1999 to 2001. 

The measurements were performed in the t:E1 channel 
at the Paul Scherrer Institute (PSI), Villigen, Switzer- 
land. The 7r+ beam with p ~ 113 MeV/c, Ap/p < 1.3 %, 
and 24 mr horizontal and vertical divergence, had an av- 
erage intensity of 6.8 • 10^7r+/s, and produced a ax,y — 
9 mm beam spot on target. A total of 2.2 • 10^^ 7r+ stops 
were recorded during the running period. 

The beam particles were first registered in a 3 mm 
thick beam counter (BC) placed ~395cm upstream of 
the detector center. The pions were slowed in a 30 mm 
thick active degrader (AD) and stopped in a 50 mm long 
segmented active target detector (AT) positioned in the 
center of the PIBETA detector. The BC, AD and AT 
detectors are all made of plastic scintillation material. 

The and /i"*" beam contaminations determined by 
the time-of-flight method were small, 0.4% and 0.2%, 
respectively. The background from hadronic interactions 
was efficiently suppressed offiine by cuts on the tim- 
ing differences between the beam detectors and the Csl 
calorimeter. 

The heart of the PIBETA detector is a spherical, 240- 
element, pure Csl electromagnetic calorimeter, 12 radia- 
tion lengths thick. The calorimeter achieves an rms en- 
ergy resolution of (Je/E = 4.0 % for 7r° mass reconstruc- 
tion at rest, sub-ns timing resolution, and an angular 
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resolution of 2.0° for minimum ionizing particles (MIPs). 
Details of the detector design, construction, performance, 
along with other relevant information regarding our ex- 
perimental method, are presented at the experiment Web 
site [lil, and in Ref. [Tg. 

Two sets of fast analog triggers accepted nearly all non- 
prompt pion decay events and a substantial fraction of 
muon decays, for energy depositions in the calorimeter 
exceeding a high threshold of HT ~ 51.7 MeV, or a low 
threshold of LT ~ 4.5 MeV. A pion stopping in the tar- 
get initiates a 180 ns long pion gate ttG, timed so as to 
sample pile-up events up to 30 ns prior to the 7r+ stop 
pulse. Event triggers were generated on the basis of a 
coincidence of the ttG gate, a beam anti-prompt veto, 
and a shower signal in the calorimeter. 
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FIG. 1: The tt"*" — > e^wy class A events (the top panel) taken 
with two-arm calorimeter trigger are virtually background- 
free (save for the 7r/3 contamination). The signal-to- 
background ratio for class B and class C events recorded with 
one-arm trigger is 3.8:1 and 7.6:1, respectively. 



Charged particles emanating from the target were 
tracked by two cylindrical multi-wire proportional cham- 
bers ^3 9.nd a 20-piece thin plastic scintillator veto ho- 
doscope (PV). The PV signals provided discrimination 
between the background protons and MIPs. MIP detec- 
tion efficiencies were continuously monitored. The av- 
erage efficiencies of the inner and outer MWPCs were 
eCi = 93.7% and = 97.9%, respectively. The average 
efficiency of the PV hodoscope was epv — 98.9%. The 
overall inefhciency for distinguishing neutral and charged 
particles was therefore 1.5 • 10~^. 

TT evj candidate events were selected from the two- 
arm HT and prescaled one-arm HT data sets by requiring 
one neutral shower in the calorimeter in coincidence with 
a positron track. In cases where there was more than one 
possible coupling (i.e., more than one neutral or positron 
track), the pair most nearly coincident in time was cho- 
sen. Our data set covers three kinematic regions: 

A: Ell\E^''^ > 51.7 MeV (two-arm HT), 

B: Elf > 20.0 MeV, E^^^ > 55.6 MeV (one-arm HT), 

C: Elf > 55.6 MeV, E^^^ > 20.0 MeV (one-arm HT), 

where the superscript "cal" refers to values measured in 
the Csl calorimeter. In all three regions the relative angle 
e'^f > 40.0°. 

The signal-to-background ratios S/B are shown in 
Fig. n where the time differences between the radiative 
photon and the positron tracks are histogrammed. The 
high energy pairs of class A show virtually no background 
{S/B > 300), while the class B (C) events with lower en- 
ergy positrons (photons) have S/B — 3.8 (7.6). 

We define the coincidence time window by Ati = |tg+ — 
^7! <5ns. The accidental background was sampled in 
two sidebands Ato, —10 ns < —t^ < —5 ns and 5 ns < 
te+ — ^7 < 10 ns. 

The accidental background is dominated by positrons 
from TT ev or the ir-fi-e decay chain, accompanied by 
an unrelated accidental neutral shower. This background 
was removed by subtracting histograms of observables 
projected using the out-of-time cut A<o from the in-time 
histograms projected via the Atj cut. 

In addition to shower energies, we also measured the 
directions of the positron and photon, thus overdetermin- 
ing the final three-body state. The positron direction was 
fixed using MWPC hits, while the photon direction was 
reconstructed from the pattern of hits in the calorimeter. 
Events with kinematics incompatible with the tt — s- evj 
decay were rejected in the final data sample. 

The time distributions of the tt — > 61^7 events with 
respect to the 7r~^ stop time are plotted in Fig. El The 
fitted exponential curves are consistent with the tt"*" life- 
time and demonstrate the purity of the final data set. 

Non- accidental background sources are (1) Tr(3 events 
for which one 7r° decay photon converts in the target, 
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producing a charged track in the detector and (2) two- 
clump showers originating from a single tt ev positron 
when a secondary shower photon or positron interacts 
in the calorimeter far enough from the primary hit to 
appear as a separate clump ("split-clump" events). 

Starting from the measured net yield of the 7r/3 decay 
events, we have used the Monte Carlo (MC) simulation of 
7 conversions in the inner detector to determine that class 
(1) background events contribute 12.9% of the signal in 
the kinematic region A and 4.3% (3.8 %) in region B (C). 
These 7r/3 contaminations were subtracted in the calcula- 
tion of TT — !■ eu^ yields. For the clump energy threshold 
E^f > 20.0 MeV used in the analysis the "split-clump" 
background (2) can be neglected. 

We have normalized the radiative pion decay rate to 
the total branching ratio of the n ev (7re2) events 
which were recorded in parallel using the prescaled one- 
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FIG. 2: tt'^ — > e"'"z^7 event timing spectra after subtraction of 
accidental background events. Predicted Monte Carlo decay 
functions are shown as full lines. 
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FIG. 3: Background-subtracted vr"'" — > e'^i/ energy spectrum 
taken with the one-arm trigger. The GEANT-simulated detector 
response is represented by the full line histogram. 



arm calorimeter trigger. The energy spectrum of the 
7r+ e^^v positrons is shown in Fig. |31 This procedure 
assures that most factors in the normalization, includ- 
ing the total number of decaying tt'^ 's and the combined 
tracking efficiency of e''"'s largely cancel out. 



The absolute branching ratio for the tt^ 
can be calculated from the expression: 



R 



7re27 



e'^v^ decay 



(3) 



where iV7re27 is the number of the detected tt — > ev^ 
events, p^e27 is the corresponding prescaling factor (if 
any), Ntj+ is the number of the decaying tt+s, = 
exp{—t/T^+)dt is the tt"*" gate fraction, A^e27 is the 
detector acceptance incorporating the appropriate cuts, 
and Ti is the detector live time fraction. An analogous 
expression can be written for the total i?7re2 branching 
ratio. Combining the two, we obtain 



R 



_ jy -^71627 P7re27 ^7re2 £-!ve2 /.n 
7re27 — ^TTe2 — Z , [^) 



where the e's denote the properly weighted products of 
epv ■ eCi • for the two data sets. The detector ac- 
ceptance ratio takes on values between 0.057 and 0.209, 
depending on the kinematic region. 

The experimental acceptances depend on both the de- 
tector response and the decay amplitudes and are calcu- 
lated in a Monte Carlo (MC) simulation. The GEANT [2^ 
simulation of the PIBETA detector response included: (i) 
the detailed geometry of the active detectors and the pas- 
sive support material, (ii) the measured detector energy 
and timing responses, (iii) event generators for tt and /j, 
decays including measured accidental pile-up rates, and 
(iv) the photo-absorption reactions in the Csl calorimeter 
incorporated in the GEANT code. 

In the Standard Model the differential rate of the 7r+ 
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6+1^7 decay can be written in the form [J]: 



-r^c2i IB {x,y) 

TT I 



Fvml 



dxdy 27r "I ' \2f.,^me 
X [ (1 + SD+ (x, y) + {\- i f SD- {x, y) ] (5) 

+ {^J^) [(1 + 7) ^it ?/) + (!- 7) ^i^t 1/)] } > 

where the /i?, SD^, and S'j^j (IB-SD interference) terms 
depend on the kinematic variables x = 2_E^/m^ and y = 
2Ee/mTr; and are the physical ("thrown") energies. 

We took the PDG value r(7re2)/rtotai = 1.230(4)- 10"'' 
1^ and calculated the experimental branc hing ratios us- 
ing the MINUIT least chi-square program [221 ■ We have 
added integral radiative corrections of —1.0% in region 
A, —1.4% in B, and —3.3% in C, to the theoretical 
^^the's |2^ . The minimization program simultaneously 
fits two-dimensional (i?g+\ E^'^^) distributions in all three 
phase space regions, constraining the integrals to the ex- 
perimental branching ratio values (i?exp)- 

The acceptance ATre2'y is recalculated in every iteration 
step of our analysis with the cuts applied to the physical 
("thrown") energies, E, following cuts applied to mea- 
sured particle energies and angles, E'^'^^, 0^^"^. Hence, our 
experimental branching ratios can be compared directly 
with theoretical absolute decay rates. 

The statistical uncertainties of the experimental yields 
are 0.6%, 1.7% and 1.5% for the regions A, B, and C, 
respectively. Systematic uncertainties of 1.8% for region 
A, dominated by 7r/3 background subtraction, 2.3% and 
3.1% for regions B and C, respectively, both dominated 
by acceptance ratio uncertainties, were added in quadra- 
ture. 

The dependence of the region-A experimental and the- 
oretical branching ratio on the value of 7 is shown in 
Fig. E (top), indicating two solutions. The positive 7 so- 
lution is preferred by ^ 50:1 once data from regions B 
and C are included in the analysis (bottom plot). We 
compare the experimental and theoretical branching ra- 
tios for the three phase space regions in Table We 
note that due to the large statistical and systematic un- 
certainties present in all older experiments, our values 
are consistent with previously published measurements. 
The best CVC fit to our data yields 7 = 0.443 ±0.015, or 
Fa = 0.0115(4) with Fy = 0.0259. This result represents 
a fourfold improvement in precision over the previous 
world average, and is consistent with chiral Lagrangian 
calculations 13, H, . 

In summary, our experimental 7r+ — s- e^z^7 branching 
ratios and energy distributions in kinematic regions A 
and C are compatible with the {V — A) interaction. The 
sizable 19 % shortfall of the measured branching ratio 
compared to the theoretical one in region B dominates 
the total x^i f^nd is disconcerting. Thus, in a fit restricted 
to region A data only, we obtain 7 = 0.480 ± 0.016, the 
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FIG. 4: Top: tv'^ e^vy branching ratio values as a func- 
tion of 7 = Fa/Fv- The theoretical parabola follows from 
the V — A model, Eq. Q. The experimental values reflect 
fits to region A data only. Bottom: minimum values of 
simultaneous fits to the entire data set (regions A, B, C). 



same as when analyzing the combined data from regions 
A and C. Significantly, all previous studies save one 
(which, too, found an anomaly [HI), have analyzed only 
data with kinematics compatible with our region A. 

Taking the {V~A) interaction and the CVC hypothesis 
as valid, the deficit could be caused by a specific detec- 
tion or analysis inefficiency in our experiment, appearing 
dominantly in region B. Detailed cross-checks, including 
absolute total tt — > ev, 7r+ — > tt^c^i^, /i — *■ eiyv, as well 
as partial ^ — > ei^Dj branching ratio evaluations from 
our data have thus far categorically excluded such inef- 
ficiencies [2^. Alternatively, the deficit could be caused 
by an inadequacy of the present V — A description of 
the radiative pion decay, i.e., Eqs. (Q), l|2Jl and (O along 
with the radiative corrections, or by an anomalous, non- 
(V-A) interaction HI IH HI . We wiU focus on the 
region B discrepancy and the latter possibility in a sep- 
arate forthcoming Letter, noting that our results clearly 
call for further theoretical and experimental work. 
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M. V. Chizhov and A. A. Poblaguev for fruitful dis- 
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TABLE I: Best-fit n —f evy branching ratios obtained with 
Fv = 0.0259 (fixed) and Fv = 0.0115(4) (fit); xVd-o.f. = 
25.4. Measured (i?cxp) and theoretical (7?the) branching ra- 
tios are shown for the three indicated phase space regions. 
Radiative corrections are included in the calculations. 
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